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I)  Introduction  AFOSR-Tfc*  8  7  -  0  74  4 

The  purpose  of  our  work  has  been  to  develop  advanced  laser 
spectroscopy  methods  to  diagnose  partially  ionized  plasmas.  We 
have  focused  on  methods  that  are  based  on  observing  the  Doppler 
shift  in  ionic  spectra  due  to  the  presence  of  an  ion  drift 
velocity.  Two  particular  methods  we  are  working  with  are 
Velocity  Modulated  Laser  Spectroscopy  (VMLS)  and  Two  Beam  Doppler 
Shift  Laser  Spectroscopy  (TBDSLS). 

The  scientific  goal  of  our  work  is  to  increase  understanding 
of  the  role  of  flow  non-uniformities  and  plasma/wall  interactions 
in  plasma  devices  by  making  in-situ  measurements  of  electric 
field  strength,  ion  mobilities,  concentrations  and  temperatures 
in  a  non-intrusi ve  fashion  that  allows  point,  one,  and  two 
dimensional  imaging. 

The  scientific  approach  i3  to  use  conventional  laser 
spectroscopic  methods  such  as  Rayleigh  scattering,  Raman 
scattering,  or  fluorescence,  to  probe  ion  absorption  line 
profiles.  If  there  is  an  electric  field  present,  the  ions  will 
experience  a  net  force  and  undergo  drift,  resulting  in  a  shift 
in  the  position  of  the  line  profile.  If  the  ion  mobility  is 
known,  then  the  electric  field  component  along  the  probe 
direction  can  be  calculated.  If  the  electric  field  driving  the 
plasma  is  modulated,  one  will  observe  an  oscillating  shift  in  the 
line  profile  that  arises  because  of  the  oscillating  force  imposed 
on  the  ions.  The  shift  may  be  related  to  the  ion  mobility,  thus 
conduct  i  v i ty . 

Temperature  and  concentration  may  be  recovered  by 
conventional  laser  spectroscopic  means.  The  methods  are  species 


x 


* 


¥ 


and  state  selective,  allowing  one  to  make  measurements  on  more 
than  one  species  and  to  study  the  effect  of  internal  mode 
nonequilibrium. 

The  merit  of  the  methods  lies  in  their  ability  to  provide 
simultaneous  measurements  of  important  parameters  in  plasmas. 

The  methods  are  well  suited  to  multi-dimensional  imaging.  One 
may  use  an  array  detector  to  image  lines  and  planes  in  addition 
to  the  more  conventional  point  configuration. 

There  has  been  a  far  amount  of  work  done  by  the  physical 
chemistry  community  on  VMLS.  VMLS  applied  to  plasmas  was 
initially  exploited  by  Woods  et  al  who  observed  Doppler  shifts  in 
microwave  transitions  of  HCO+  [1,2]  and  then  used  the  shifts  to 
infer  ion  mobility  [3].  Microwave  measurement  of  HCO+,  NO+,  C0  + 
and  HNN*  were  performed  by  DeLucia  et  al  [4-7]. 

VMLS  was  extended  to  the  infrared  region  by  two  groups,  that 
of  Professor  Richard  Saykally  of  the  Chemistry  Department  at  UC 
Berkeley  [8-14]  and  that  of  Professor  T.  Oka  of  the  Chemistry 
Department  at  the  University  of  Chicago  [15-19].  Say’  a  ly's  group 
has  focused  on  the  spectroscopy  of  ions  in  plasmas  while  Oka  has 
focused  on  measurements  of  ion  mobility. 

The  Chicago  group  has  observed  Ar H + ,  NH4,  HCNH+  and  HC0P  +  . 

The  Berkeley  group  has  made  observations  of  H^O*.  NHjJ,  H2F*. 

HCO\  HN N *  ,  and  H3  • 

The  Berkeley  group  has  extended  the  method  to  the  visible 
region  [20-21].  They  have  observed  both  C0+  and  N2  in  CO  or  N2 
seeded  He  plasmas. 

The  probe  method  in  all  the  studies  so  far  is  absorption 
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spec toacopy .  The  infrared  work  waa  done  with  a  diode  laser 
(Chicago  group)  and  a  Krypton  pumped  color  center  laser  (Berkeley 
group.)  The  visible  work  was  done  with  a  single  mode  dye  laser. 
One  of  the  unique  aspects  of  our  work  is  that  we  are  exploring  a 
number  of  different  probe  methods  so  as  to  be  able  to  extend  the 
method  into  the  ultraviolet  and  to  molecules  that  do  not  absorb 
at  appropriate  frequencies. 

In  the  following  we  describe  the  theoretical  basis  of 
doppler  shift  methods,  outline  the  nature  of  the  two  methods  we 
are  pursuing,  and  describe  experimental  results  to  date.  The 
work  is  then  summarized  and  appropriate  conclusions  drawn. 


1 
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II)  Theoretical  Basis 


a)  Ion  Response  to  a  Modulated  Electric  Field 

The  basis  of  Doppler  shift  methods  is  Coulomb's  Law  which 
states  that  the  force  felt  by  a  charged  particle  is  directly 
propotional  to  the  imposed  electric  field 

F  -  q  E 

(1  ) 

where  E  is  the  electric  field,  q  the  charge  of  the  particle  and  F 
the  force.  In  a  collisional  plasma,  one  may  write  a  momentum 
equation  for  the  ions  that  is  of  the  form  (in  the  absence  of  a 
magnetic  field) 

d*/d  t  -  (q/m)E  -  voT 

(2) 

where  m  is  the  ion  mass,  and  vc  the  collision  frequency.  IT 
there  is  a  DC  electric  field  present,  then  the  ion  drirt  velocity 
will  be  (after  an  initial  transient) 


*  •  Mi  E 


where 


-  q/m vc 


is  the  ion  mobility.  Likewise,  if  the  electric  field  is 
oscillating  in  form 


(3) 


E  «  E  exp(iwt) 

(  U  ) 

where  w  is  the  angular  frequency,  ,  then  the  velocity  will 
be  modulated  (after  an  initial  transient)  in  the  form 
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T  "  ( U^/(l +i(u/vc)))  B0  exp(lwt) 

(5) 

Because  of  the  Doppler  effect,  the  observed  resonant 
absorption  or  emission  frequency  of  a  molecule  with  a  non-zero 
drift  velocity  is  "hifted  according  to  the  relation 


dv  -  ^0/c)  v 

(6) 

where  vQ  i3  the  line  center  frequency,  c  the  speed  of  light,  and 
v  the  component  of  velocity  along  the  line  of  sight.  Thus,  for  a 
DC  field,  the  shift  will  be 


dv/v  -  »iE/c 

(7) 

and  the  depth  of  modulation  in  resonant  frequency  due  to  an 
oscillating  field  will  be 


dv/vQ  -  ( ( u./c ) /( 1  +  i ( w/ ve))  E0  exp(iwt) 

(8) 

where  E  and  E0  are  the  component  of  the  electric  field  along  the 
line  of  sight  in  the  DC  and  AC  cases  respectively.  The  magnitude 
of  the  AC  shift  will  be 


dv/vQ  -  (p^/c)  (  1 ♦ ( u / vc ) 2 ) ~  1  / 2  E0 

(9) 

and  the  phase  of  the  shift  will  be 

$  ■  tan*' (-u/Vp) 

(10) 

Note  that  in  the  limit  of  small  w/vc>  the  modulated  result 


becomes  the  same  as  the  DC  result.  When  u>>>vC(  however,  the 
modulation  depth  goes  to  zero  and  is  independent  of  the  collision 
frequency.  Since  the  experimentalist  controls  w,  however,  this 
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is  not  a  serious  problem.  If  the  collision  frequency  were  j 

extremely  small  the  plasma  would  be  in  a  collisionless  domain.  i 

i 

In  practice,  ions  in  a  hot  plasma  will  have  a  thermal 
velocity  distribution  that  will  result  in  a  fairly  broad 
absorption  line  profile.  The  modulated  shift  will  be  easier  to 
see  if  it  is  large  with  respect  to  the  profile.  If  we  take  the 
line  width  to  be  given  approximately  by  the  Doppler  profile  we 
can  estimate  the  ease  with  which  a  shift  can  be  seen.  The 
equilibrium  Doppler  width  is 


dvD/v0  -  (Zln2/c)  vt 


(11) 


where  vt  is  the  mean  thermal  velocity.  Since  the  ion  collision 
frequency  may  be  written 


vc  -  nTQvt 


(12) 


where  n?  is  the  total  number  density  and  Q  the  collision  cross 
section,  one  may  write  (if  u  <<  vc  0r  for  the  DC  case) 


dv/dvjj  .  q/ ( (/ln2 )n jmQvf )  E 
Invoking  the  perfect  gas  law  and  noting  that 


(13) 


vt  -  (8kT/*m)1/2 

(1  M 

where  k  is  Boltzmann's  constant  and  T  the  ion  temperature,  we 
obtain 


dv/dvD  -  ( xq/ 8Q/l n  2 )  (E/P) 


where  P  is  the  gas  pressure. 


(15) 
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The  coefficient  (irq/8Q/ln2)  is  of  order  unity  in  typical 
plasaas  for  E  in  volts  per  meter  and  pressure  in  Pascals.  Thus 
one  can  obtain  significant  modulation  depths  or  DC  shifts  at 
moderate  values  of  E. 

b)  The  Effect  of  a  Magnetic  Field 

The  effect  of  a  steady  magnetic  field  is  to  cause  the 
conductivity  to  become  a  tensor  quantity.  In  the  direction  of 
the  magnetic  field,  however,  the  conductivity  is  the  same  as  that 
given  above.  Thus  by  a  suitable  experimental  arrangement  one  may 
still  utilize  the  above  equations.  Furthermore,  if  the  collision 
frequency  is  large  compared  to  the  plasma  gyro  frequency,  then 
the  tensor  conductivity  reduces  to  its  zero  magnetic  field  scalar 
limit. 

The  plasma  gyro  frequency  is 


<*>  r  »  q  B  /  m 

(16) 

For  singly  ionized  argon,  the  value  of  q/m  is  about  5.32E6  C/kg, 
so  that  for  a  magnetic  field  of  one  Tesla,  the  gyro  frequency 
would  be  of  the  order  of  a  megahertz.  Under  collisional 
conditions  in  most  plasmas,  the  collision  frequency  will  be 
higher  than  this.  If  not,  the  full  equations  must  be  used, 
c)  Spectroscopic  Methods 

The  choice  of  a  spectroscopic  method  depends  on  the 
spectroscopic  properties  of  the  ionic  species  to  be  probed  and  on 
the  thermodynamic  state  of  the  plasma. 

The  candidate  methods  include  single  and  multi-photon  laser 
induced  fluorescence,  Rayleigh  scattering  and  electronic  Raman 
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scattering  (both  coherent  and  incoherent.)  If  molecular  ions  are 
of  interest,  then  rotational  and  vibrational  Raman  scattering  are 
possible  candidates. 

At  the  present  time  we  must  divide  ionic  species  into  two 
spectroscopic  classes  according  to  whether  their  lowest  lying 
electronic  states  can  or  cannot  be  directly  excited  by  available 
laser  systems.  Hydrogen  and  argon  are  two  common  propelants 
whose  ionic  ground  states  are  extremely  difficult  to  access 
directly.  Barium  is  an  element  commonly  added  as  easily 
ionizable  seed  with  an  easily  accessed  ionic  ground  state. 

The  thermodynamic  state  of  the  plasma  is  important  in 
determining  the  possible  modes  of  operation.  For  a  plasma  that 
i s  fully  ionized,  one  may  utilize  non-resonant  methods  such  as 
Rayleigh  scattering  because  the  plasma  is  almost  pure  ions  and 
electrons.  If  the  degree  of  ionization  is  much  less  than  unity 
or  there  is  more  than  one  ionic  specie  present,  one  must  be  able 
to  distinguish  between  species,  requiring  a  resonant  method  such 
as  laser  induced  f lourescence. 

d)  Signal  Behavior  for  VMLS 

Consider  the  case  in  which  we  chose  a  spectroscopic  method 
(laser  induced  flourescence  or  Rayleigh  scattering,  for  example 
see  Figure  1)  that  probes  the  Doppler  profile.  That  is  for  which 
the  signal  is  proportional  to  the  Doppler  line  shape  function. 

If  the  laser  source  has  a  line  width  much  narrower  than  the 
Doppler  width  then  the  signal  will  be  directly  proportional  to 
the  Doppler  line  shape  function 


Sig  «  (2/ln2/AvD/,r)  exp-((2/ln2/AvD)( v-vQ))2 


(17) 


The  effect  of  velocity  modulation  on  the  Doppler  profile  is 
felt  as  a  modulation  of  v0.  By  differentiating  the  signal  twice, 
one  obtains  the  result  that  the  signal  is  most  sensitive  to 
changes  in  vQ  when  vl,  =  vQ  ±  O.M25  Avp  and  the  fractional 
sensitivity  is 

dSig/Sig  =  2.35  dvn/avn 

(18) 

Using  the  expression  for  the  modulation  depth  previously  derived, 
we  may  recast  this  expression  into  a  detectabilty  limit  for  the 
detection  of  a  given  modulation  voltage 

EoDet  “  1*08  (Q/Q)  P  (dSig/Sig)Det 

If  one  assumes  a  moderate  one  percent  for  the  signal 
detectability  limit,  and  a  typical  ten  square  Angstroms  for  the 
collision  cross  section,  then  one  obtains 

EoDet(v/ra>  -  0.00675  P(Pa) 

(20) 

At  one  atmosphere,  the  detectable  electric  field  is  only  6.75 
V/cm,  a  field  readily  achievable  in  the  laboratory.  By  utilizing 
phase  sensitive  detection,  one  can  achieve  significantly  lower 
detectability  limits. 

d)  Signal  Behavior  for  TBDSLS 

In  principle,  any  method  which  probes  the  Doppler  profile 
with  sufficient  resolution  can  be  used  to  determine  a  DC  line 
shift.  If  o  ly  one  probe  is  used,  however,  there  may  be  a 
problem  of  absolute  frequency  calibration.  By  using  two  beams  in 


opposite  directions,  an  in-situ  calibration  is  achieved. 

Consider  the  arrangement  illustrated  in  Figure  2.  Assume 
that  there  is  a  drift  velocity  in  the  positive  x  direction.  If 
the  laser  frequency  is  scanned  through  the  line,  then  the 
rightward  running  beam  will  probe  the  profile  labeled  R  in  figure 
2,  while  the  leftward  running  wave  will  probe  the  profile  labeled 
L. 

As  the  laser  frequency  is  scanned  through  the  line  (or  lines 
as  it  were,)  from  v  <  vQ  to  v  >  v0,  signal  will  first  appear  from 
the  right  running  beam.  As  the  frequency  is  increased,  the 
leftward  running  beam  will  begin  to  contribute.  If  the  two 
signals  are  independently  observed,  they  will  look  like  that  of 
Figure  2.  (By  chopping  the  beams  at  a  rate  fast  compared  to  the 
frequency  scan  rate,  one  may  use  the  same  imaging  optics  to 
detect  both  signals.)  The  peak  separation  is  just  twice  the 
Doppler  shift  and  the  crossing  point  of  the  two  signals  is  the 
unshifted  line  center. 

For  the  method  to  be  useful  one  must  also  simultaneously 
determine  the  ion  collision  frequency.  If  one  assumes  to  first 
order  that  the  collision  cross  section  is  constant,  then 


vc  oc  p//t 

(21) 

Thus,  one  may  obtain  relative  profiles  without  knowing  the  cross 
section,  or  absolute  profiles  if  calibration  is  possible  or  the 
cross  section  is  known.  The  temperature  can  be  obtained  by 
fitting  the  line  shape  of  one  of  the  lines  to  the  appropriate 


Voigt  or  Doppler  profile,  or  by  other  means  when  the  line  profile 


Is  insensitive  to  temperature.  The  pressure  can  be  measured  by 
conventional  means. 

The  detectability  limit  for  the  method  is  determined  by  the 
resolution  with  which  the  distance  between  the  two  peaks  can  be 
resolved.  If  the  precision  with  which  the  line  shape  signal  can 
be  measured  is  As,  then  for  Doppler  broadening  the  precision  with 
which  the  separation  in  peaks,  and  thus  electric  field,  can  be 
resolved  is 


AE/E  *  ( 1  //21n2 )  (A\>q  Avg)  (As/s)^^ 


If  we  define  the  detectability  limit  as  that  value  of  E  for 
which  a  precision  of  10%  is  achieved  when  the  spectroscopic 


signal  is  detected  with  a  precision  of  1  % ,  then  we  can  write 
approximate! 


Edet(V/cm)  -  0.04  PUarJ/TOO1  /2 

(23) 

(Assuming  that  Q  «  10A2,  ra  is  the  mass  of  a  proton,  and  the  line 
center  is  at  5000A.)  Thus  at  a  pressure  of  0.01  Bar  and  5000K, 
the  detectability  limit  is  about  1  V/cm. 

It  should  be  noted  that  least  squares  fitting  the  line  shape 
to  the  appropriate  profile  as  suggested  above,  would 
significantly  increase  the  precision  with  which  the  line  center 
is  detected  and  thus  reduce  the  detectability  limit  for  the 
electric  field.  Also,  if  a  suitable  transition  is  found, 
saturation  spectroscopy  [22]  could  be  used  to  locate  the  line 
peak  with  high  precision,  although  the  in  situ  calibration  for 
the  zero  shift  position  would  be  lost. 


e)  Effect  of  Plasma  Conditions 

Plasmas  of  Interest  in  propulsion  applications  span  a  wide 
range  of  state  conditions.  At  the  inlet  of  a  MPD  thruster,  for 
example,  the  pressure  is  quite  high  and  the  gas,  as  yet  un¬ 
ionized,  is  collision  dominated.  As  the  gas  becomes  ionized  and 
the  resulting  plasma  is  accelerated  and  expands  through  the 
exhaust  nozzle,  the  pressure  drops  dramatically  until  collisions 
become  unimportant.  The  temperature  may  start  out  very  low,  room 
temperature  in  the  laboratory,  but  can  rise  to  values  as  high  as 
fifty  or  one-hundred  thousand  degrees  Kelvin. 

Pressure  and  temperature  changes  affect  the  diagnostic 
mainly  through  changes  in  the  collision  frequency  and  in  the 
ratio  of  modulation  depth  to  the  Doppler  width.  As  long  as  the 
former  is  large  compared  to  the  gyro  frequency  and  the  latter 
large  enough  to  detect,  the  diagnostic  will  operate. 

Figure  3  is  a  parametric  plot  of  the  collision  frequency, 
the  ratio  of  modulation  to  Doppler  width,  and  the  mean  free  path 
as  a  function  of  pressure  and  temperature.  A  collision  cross 
section  of  ten  square  Angstroms  and  an  electric  field  of  one 
hundred  volts  per  centimeter  were  used  to  calculate  the 
parameters.  As  can  be  seen,  as  the  pressure  drops,  the 
modulation  depth  increases  at  a  much  faster  rate  than  the 
collsion  frequency  decreases.  Thus  the  flow  becomes 
collisionless  long  before  the  measurement  fails. 

Another  consideration  is  pressure  or  Stark  broadening  at 
high  electric  field  strengths  .  The  methods  will  still  work  even 
if  pressure  or  Stark  broadening  [23]  make  an  important 
contribution  to  the  line  width,  although  they  are  best  suited  to 


the  Doppler  dominated  and  low  ion  and  electron  concentration 
limits.  The  detectability  limit  will  merely  be  increased  because 
it  is  determined  by  the  peak  detection  resolution.  The  broader 
peak  of  a  pressure  or  Stark  broadend  line  is  harder  to  resolve. 

The  Stark  effect  will  cause  shifts  in  the  line  position  due 
to  both  the  imposed  electric  field  and  the  field  of  electrons  and 
other  ions.  The  DC  Stark  effect  can  be  accounted  for  in  a 
straightforward  fashion,  and  at  some  plasma  conditions  will 
dominate  the  shift. 

The  presence  of  a  Stark  shift  due  to  electrons  and  ions  will 
be  more  of  a  limiting  factor.  Unless  one  can  calibrate  out  the 
shift,  it  will  appear  as  a  systematic  error.  Of  course  if  the 
line  width  is  dominated  by  Stark  broadening,  then  the  Stark  width 
can  be  measured  and  the  contribution  of  the  Stark  effect  to  the 
measured  shift  calculated. 


Ill)  Experimental  Apparatus 


Our  system  is  illustrated  in  Figure  4.  We  have  an  argon  ion 
laser  (Spectra  Physics  171-19)  pumped  ring  dye  laser  (Coherent 
699-21).  The  ring  laser  is  fully  stabilized  and  has  a  line  width 
of  less  that  one  megahertz.  The  output  of  the  laser  is  directed 
through  a  Fabry-Perot  interferometer  (Coherent  Model  251  Spectrum 
Analyzer)  as  a  means  of  measuring  the  relative  wavelength 
precisely.  A  small  portion  of  the  beam  is  also  split  off  and 
sent  to  a  wavemeter  (Burleigh  Model  WA  10)  for  absolute 
wavelength  calibration. 

The  beam  is  then  directed  into  the  test  section  (Figure  5) 
in  either  the  VMLS  or  TBDSLS  configuration.  In  either  case,  the 
beam  (or  beams)  pass  through  small  holes  in  the  electrodes.  A 
power  supply  can  deliver  several  hundred  volts  DC  modulated  at  up 
to  100  VAC. 

The  signal  is  detected  by  observing  the  fluorescence  at 
ninety  degrees  (Figure  4).  The  signal  is  collected  with  F/9 
optics  and  passed  through  a  1/4  meter  Jarrell-Ash  monochromator. 

A  RCA  1 P 2 8  photomultiplier  is  used  as  the  detector. 

The  plasma  used  in  the  proof  of  principle  experiments  is  a 
flame  generated,  partially  ionized  plasma  in  which  the  dominate 
ion  is  barium.  Barium  was  chosen  as  the  test  ion  because  of  the 
low  ionization  potential  of  the  neutral  and  because  it  has 
a  convenient  structure  with  readily  accessible  absorption  lines. 
We  pump  at  455  nm  and  detect  c o  1 1  i  s  i  o n a  1  1  y  induced  fluorescence 
at  493.5  nm  to  avoid  scattering  interference. 

The  burner  is  a  capillary,  diffusion  flame  burner  designed 


by  Krupa,  at.  al.  [24].  It  allows  a  wide  range  of  operating 
conditions  without  the  normal  difficulty  of  using  acetelyne  and 
oxygen  as  reactants.  The  burner  Is  mounted  In  a  low  pressure 
vessel  and  can  be  operated  down  to  about  20  Torr.  Barium 
tetrachloride  Is  seeded  Into  the  flame  using  a  Perkin-Elmer 
aspirator  modified  to  operate  lower  pressures. 


IV)  Progress 


Our  grsnt  initiated  in  December,  1985.  During  the  first 
year  of  the  grant  we  have  developed  the  theory  for  the  methods 
and  started  experiments  to  demonstrate  Velocity  Modulated  Laser 
Spectroscopy  (VMLS)  and  Two  Beam  Doppler  Shift  Laser  Spectroscopy 
(TBDSLS ) .  Calculations  indicate  that  the  methods  will  lead 
to  excellent  detectability  for  measurement  of  electric  field  (or 
ion  mobility),  species  concentration,  and  temperature. 

We  also  observed  fluorescence  from  the  barium  ion  and 
measured  the  *»55  no  line  width  at  atmospheric  pressure.  At 
atmospheric  pressure  the  linewidth  is  dominated  by  pressure 
broadening.  Figure  5  shows  a  typical  line  shape  profile.  The 
series  of  sharp  peaks  is  the  Fabry-Perot  signal  superimposed  to 
provide  a  measure  of  the  frequency.  The  peaks  are  1.5  GHz  apart, 
thus  the  line  is  about  A. 5  GHz  wide.  We  have  revised  the 
detectability  calculations  to  account  for  the  pressure 
broadening.  We  started  the  new  year  conducting  experiments  at 


low  pressures. 
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